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ABSTRACT 


Motivated by recent visits from interstellar comets, along with continuing discoveries of minor bodies 
in orbit of the Sun, this paper studies the capture of objects on initially hyperbolic orbits by our solar 
system. Using an ensemble of ~ 500 million numerical experiments, this work generalizes previous 
treatments by calculating the capture cross section as a function of asymptotic speed. The resulting 
velocity-dependent cross section can then be convolved with any distribution of relative speeds to 
determine the capture rate for incoming bodies. This convolution is carried out for the usual Maxwellian 
distribution, as well as the velocity distribution expected for rocky debris ejected from planetary 
systems. We also construct an analytic description of the capture process that provides an explanation 
for the functional form of the capture cross section in both the high velocity and low velocity limits. 


Keywords: Solar system (1528), Dynamical evolution (421), Small solar system bodies (1469), Kuiper 
belt (893), Oort cloud (1157) 


1. INTRODUCTION 


The past few years have witnessed the detection of two interstellar bodies passing through the solar system on 
hyperbolic orbits. The discoveries of the irregular body ‘Oumuamua (Meech et al. 2017) and the comet Borisov 
(Jewitt & Luu 2019) sparked immediate interest in characterization of these objects and facilitated wide-ranging 
speculation regarding the possibility that our solar system is more broadly contaminated by minor bodies of extra- 
solar origin (e.g., Siraj & Loeb 2019; Namouni & Morais 2020). Although no current evidence indicates that any 
specific objects in the solar system are of extrinsic origin (Morbidelli et al. 2020), the question of whether or not any 
such objects reside in interplanetary or trans-Neptunian space is of considerable interest. Motivated by these issues, 
this paper reconsiders the capture of external bodies by our solar system. The calculation of the capture cross sections 
is the first step in assessing whether or not the solar system presently contains quasi-permanently trapped interstellar 
bodies. This treatment also provides constraints on the expected orbits of any such material. 

The dynamics of the outer solar system represents one of the oldest problems in theoretical astrophysics. Starting 
more than two centuries ago, classic studies include the long-term stability of the solar system (Lagrange 1776; Laplace 
1799-1825), the origin of comets (Laplace 1846), and orbital anomalies that led to the discovery of Neptune (Le Verrier 
1846; Adams 1846). Over recent decades, the outer solar system has revealed itself to be increasingly complicated, 
with the discovery of the Kuiper Belt (Luu & Jewitt 2002), dozens of dwarf planets (starting with Brown et al. 2004, or 
perhaps Tombaugh 1946), high-inclination objects (Becker et al. 2018), and aligned extreme trans-Neptunian objects 
(Sheppard & Trujillo 2016) that led to the hypothesis of a possible ninth planet (Batygin & Brown 2016; Batygin 
et al. 2019). The more recent discovery of interstellar objects (Meech et al. 2017; Jewitt & Luu 2019) adds to the 
intrigue. Both the complex orbital architecture of the solar system and the presence of interloping objects motivates 
this present study. The goal is to determine cross sections for the capture of foreign bodies by the solar system, and 
to obtain a deeper understanding of the capture process. 
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The possible capture of interstellar bodies by the solar system also has a long history. The general problem of 
interacting binaries was considered by Heggie (1975), where the subset of ‘resonant’ encounters lead to capture. 
Subsequent studies have carried out numerical explorations of the capture process specifically for our solar system, 
often considering only the Sun-Jupiter system (see, e.g., Valtonen & Innanen 1982; Valtonen 1983; Siraj & Loeb 2019). 
Additional studies consider capture for specific scenarios, including capture by compact objects (Pineault & Duquet 
1993), capture of interstellar objects from the field (Lingam & Loeb 2018; Hands & Dehnen 2020), the formation of 
wide binaries (Kouwenhoven et al. 2010), and the possible capture of Planet Nine (Li & Adams 2016; Mustill et al. 
2016). Most of these previous studies calculate the capture rate by sampling a given distribution of encounter speeds 
between the incoming body and the solar system. These studies generally use the field star velocity distribution, with 
dispersion ~ 40 km/s (Binney & Tremaine 2008) or that appropriate for the solar birth cluster (Portegies Zwart 2009; 
Adams 2010; Pfalzner 2013; Parker 2020), where the velocity dispersion is expected to be ~ 1 km/s (Lada & Lada 
2003). Notice, however, that the velocity distribution for rocks (or planets) ejected by solar systems will not generally 
have a simple Maxwellian form.! 

The objective of this paper is to extend the aforementioned previous work concerning the capture of interstellar 
bodies by the solar system. Whereas most studies determine capture rates and cross sections for a given distribution 
of velocities, this work finds the cross section ø (væ) as a function of relative velocity. The results can then be integrated 
(after the fact) for any distribution of velocities of interest. This approach is much more computationally expensive than 
previous treatments, but is made possible with current computational capabilities. Specifically, this paper reports the 
results from ~ 5 x 108 fly-by simulations. In addition, we carry out the simulations for solar systems models including 
all four giant planets. Although earlier work (Heggie 1975; Pineault & Duquet 1993) provides analytic estimates for 
the cross sections, exact forms are not available (primarily due to the lack of an analytic solution to the gravitational 
three-body problem). We revisit this issue using a different (but equivalent) set of approximations. We then compare 
the numerical and analytic results for the cross section as a function of velocity, and find good agreement. 


2. DYNAMICS OF THE ROCK CAPTURE PROCESS 


This section presents an analytic description of the rock capture process. The capture of an incoming body occurs 
through the time dependence of the gravitational potential of the solar system. In this treatment, we consider the 
incoming orbit in two regimes. In the outer regime, at large distances, the rock executes a hyperbolic orbit about 
the center of mass of the solar system. In the inner regime, at closer distances, the rock can enter into the sphere of 
influence of individual solar system members (e.g., the Sun or Jupiter), and then be described by a hyperbolic orbit 
around that body. Under favorable conditions, the deflection by the solar system body during the close encounter can 
lead to energy loss and capture in the center of mass frame. This effect is essentially the inverse of the gravitational 
slingshot mechanism by which satellites are boosted through planetary encounters. Note that by dividing the orbit 
into two regimes, we are implicitly assuming that 3-body effects are not important. 

For the sake of definiteness, we consider only one planet at a time, and work in the limit where the masses of the 
rock u, the planet m, and the star M obey the ordering 


uaeme A, (1) 


The incoming orbit of the rock is characterized by its asymptotic speed vœ and impact parameter b. For given input 
variables (v..,b), we can define the orbital elements and related physical quantities, including the specific energy and 
angular momentum, 


1 
E= SN and J = bn, (2) 
the semi-major axis and eccentricity, 
GM 
la| = —a = 2 and =e? =1+0?/a?, (3) 
and the perihelion distance 
Tp =p=a(1—e) =|al(e—1). (4) 


l As one example, the distribution of speeds for planets ejected from crowded solar systems has the approximate form dP/dv = 4u/(1 +?) 


(e.g., Moorhead & Adams 2005). 
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Note that, by convention, the semi-major axis a < 0. To fully characterize the orbit, one must also specify the 
inclination angle of the incoming trajectory. 

It is useful to define the effective cross section for hyperbolic orbits to enter the giant planet region of the solar 
system. In order for the incoming rock to experience the time-dependence of the gravitational potential, the perihelion 
rp must be smaller than the semi-major axis a, of the planet of interest. This condition implies that the impact 
parameter b is bounded from above by b? < a? + 2a,|a|, where a is the semi-major axis of the incoming orbit. The 
nominal cross section ou for orbit crossing is thus given by 


GM 
ee: , (5) 


Co 


ou lo + 2a,|a|] = 27 


where the final equality holds for essentially all incoming speeds of interest (v2, < GM/ap). The capture cross section 
will be some fraction of the fiducial cross section (5).? 


2.1. Gravitational Slingshot Mechanism for Close Encounters 


For the inner regime defined above, we consider close encounters of incoming rocks on initially hyperbolic trajectories 
with much larger target bodies (either the Sun or one of the giant planets). We can define the coordinate system so 
that the rock approaches the target body from the +4 direction and from the + direction, where the angle of the 
incoming trajectory is 0 in the center of mass frame (see the first panel in Figure 1). The rock initially has speed D 
in the inertial reference frame and the target body has speed U. In the frame of the target, the incoming rock has 
velocity 


Vitar = (—[vz = U], vy) (6) 


and the outgoing velocity has the form 
Vatar = (+[v2 — U], vy) « (7) 


This second equation assumes that the encounter is symmetric, i.e., the outgoing trajectory of the rock is the mirror 
image of the incoming trajectory. This approximation thus assumes that the larger body does not change its velocity 
(consistent with the ordering of equation [1]) and that the encounter time is short compared to the orbital period. In 
the center of mass reference frame, the incoming velocity has the form 


Viem = (—vz, vy) = (—vcos 0, vsin 0) , (8) 
whereas the outgoing velocity becomes 
Vaem = (+[vx — 2U], vy) = (+[v cos 0 — 2U], v sin 0) . (9) 
The final speed is then given by the expression 
vs = [v cos — 2U]? + v? sin? 0 = v? — 4Uucos@ + 4U? . (10) 


Note that this discussion assumes that the encounter is symmetric in the reference frame of the capturing body. This 
approximation is expected to be valid because only close encounters with the body result in capture events, and such 
close encounters will be symmetric to leading order. 

The discussion thus far has implicitly assumed that the target is moving in the —ĉ direction. In general, however, 
the target can also have a ĝ component to its velocity. Because of the geometry of the encounter, however, only the 
& component of the rock velocity changes (in this approximation). We can thus incorporate the more general case by 
interpreting the velocity U as the component of the target velocity in the —ĉ direction. With this definition of U, the 
final speed still obeys equation (10). 


2 Note that the interpretation of this fiducial cross section would be more complicated if the planetary orbit had significant eccentricity. 


Nonetheless, one can always scale the results to the expression of equation (5). 
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Figure 1. Examples of capture events. In each frame, the dotted black line denotes the rock’s initial (unbound) orbit; the 
solid black line denotes the rock’s trajectory after a perturbation; and the arrows specify the direction of the orbit. The yellow 
circle represents the Sun, the filled red circle represents Jupiter’s sphere of influence drawn at the epoch of the rock’s closest 
approach, and the red line marks Jupiter’s orbit. The top row shows a capture by Jupiter: the left panel is in the frame of the 
solar system’s barycenter, and the right panel is in Jupiter’s rest frame. Note that the rock gets well inside of Jupiter’s sphere 
of influence, but does not actually collide with the planet. The bottom row shows a capture by the Sun: the left panel is in 
the frame of the solar system’s barycenter, and the right panel is in the Sun’s rest frame. If the target (Sun or planet) has a 
component U of its velocity moving away from the incoming rock as it approaches periapsis, then the encounter causes the rock 
to lose energy in the inertial reference frame, thereby allowing the rock to potentially enter into a bound orbit. 


2.2. Solar Close Encounters 
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When the rock is far away from the Sun, it effectively orbits the center of mass of the system. When the radial distance 
of the rock becomes sufficiently small, however, its orbit is determined by the location of the Sun. We can delineate 
the boundary between these two regimes by considering the acceleration (and hence forces) in an accelerating reference 
frame. As expected, the orbit of Jupiter provides an approximate boundary between the outer problem (hyperbolic 
orbit about the center of mass of the solar system) and the inner problem (close encounter with the Sun). 

When the rock enters the sphere of influence of the Sun, its speed is given by 

2 2 , 246M 


H KE (11) 


where r < ay is the location of the rock. It will then execute a (hyperbolic) orbit around the Sun. Due to the motion 
of the Sun about the center of mass of the solar system, the post-encounter velocity will be given by 


vs = v? — 4Uv cos 0 + 4U", (12) 


where U is the component of the solar velocity in the direction of the perihelion of the orbit and @ defines the shape 


of the hyperbola. Capture of the rock requires that v3 < v? — v2,, so that we obtain the constraint 


v2, < 4Uvcos@ — AUT. (13) 
Here, the angle 0 is determined by the parameters of the original hyperbolic orbit about the Sun, so that 
1 
cos = - = saa (14) 


e (a2+b2)1/2’ 


where a is the semi-major axis and b is the impact parameter. The capture constraint thus becomes 


cn 1/2 1/2 2 
5 m (GM 2 , 2GM ja] i m 
Voo <ta (=) D Tp (a2 + 62)1/2 -O M2)’ (15) 


where we have assumed that the speed U is some fraction of the speed of the Sun in its orbit about the center of mass. 
Specifically, we define the parameter o such that 


Uza o (16) 


where m is the mass of Jupiter. Since the speed of the incoming rock v >> U for the close encounters of interest, 
we ignore the U? term. Finally, we evaluate the rock velocity at the perihelion position rp (see equation [4]), as this 
location corresponds to where the close encounter takes place. Working to consistent order, the expression for the 
capture constraint can be written in the form 


v2 < 4a (CM 1/2 T4GM]a]] P ` ga (CM 1/2 GM an 
a M \ aj Ve "gr buso ` 
The constraint can be written as a limit on the impact parameter, i.e., 
m (GM\3/? D ES 
b < 8a— | —— org 93 — Í 18 
Sch ay ) Yoo wa (; =) GE 


If one requires that the rocky body is not only captured, but is captured into an orbit with semi-major axis less 
than some maximum value amax, then the left-hand-side of equation (17) can be replaced with v2, + v2, where v2 = 
GM/dmax- Finally, note that this treatment implicitly assumes that U > 0. If the Sun is moving in the opposite 
direction, the encounter would cause the incoming rocky body to gain energy, and capture does not take place. 
Given the approximations presented above, the resulting cross section for capture can be written in the from 


2/(GM\? 1 
= Biesen a2 (=) 
o = Dao ox M F (G2. eye (19) 


This cross section is specified up to the constant a, which is expected to be of order (but less than) unity. This form 
is consistent with those derived earlier by other means (Heggie 1975; Pineault & Duquet 1993; Valtonen 1983). Notice 
that this derivation breaks down for sufficiently high incoming speeds, Ga 28 km/s, as shown in Appendix A. 
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2.3. Planetary Close Encounters 


Another channel for capture occurs through close encounters with the giant planets, most often Jupiter, which will 
be considered in this discussion. Equation (5) represents the cross section for an incoming rock to enter the sphere of 
radius ay. Only a fraction of the incoming trajectories fı = R2,,/4a%, will enter the sphere of influence of Jupiter,’ 
delineated by Rs; ze az(m/M)?/5 (Bate et al. 1971). However, not all of the orbits that enter the sphere of influence 
will pass close enough to the planet to experience significant deflection. As a result, we must estimate the smaller 
fraction fo of trajectories that allow for capture. 

As a rough approximation, significant deflection requires cos to be of order (but still less than) unity, which in 
turn implies bap ~ |a|np (equation [14]), where (anp, bap) correspond to the elements of the hyperbolic orbit around the 
planet. When the rock encounters the planet, its speed in the solar reference frame is given by equation (11) evaluated 
at r+ ay. The asymptotic speed (Voo)ap for the hyperbolic orbit about the planet depends on the planetary motion, 
but will typically be of the same order. We can thus write 

CB) = 0, + 9 =v} +02, (20) 
where ( is a dimensionless factor of order unity and where the second equality defines the velocity scale v,. The 
semi-major axis of the hyperbolic planetary encounter is given by 

Gm Gm m 


(la|)ap = ( = 


Ce EE 
ul vU2+u2 M 


(21) 


Since we require bap Ș |a|np and |alnp < Rsor, the fraction fz = Gin» /4a7,. The resulting cross section for capture due 
to planetary encounters has the form 


_ om GM Gm  \°_ yma a Sek 1 Ge 
Oo we, Taye ail 2 \M ay) (v2, +02)?’ 


co 


where we have introduced a dimensionless factor y that is expected to be of order unity. Note that this expression has 
a form similar to that of equation (19), which corresponds to the capture cross section for solar encounters. Keep in 
mind, however, that the velocity scales are different and are expected to obey the ordering vy < vz. 


2.4. Energy Distribution of Newly Bound Orbits 


Using the results from the previous section, we can write the post-encounter speed of the rock in the form 
vgl — 4Uvcos0. (23) 


The semi-major axis a, of the bound orbit is defined so that 


GM GM 
=E =2 — vs = 4Uv cos 0 — v2,. (24) 
Or r 
Let us now define a scale length bọ according to 
m (GM\?” ay ae 
bo = 100 a 25 
£ M ( ay ) Bee Se Lë ES 
With this construction, the semi-major axis of the bound orbit is given by 
GM bo agb aox 
IER b = bn bh 1-x’ (26) 


where we let ao = |a| denote the (magnitude of) the semi-major axis of the initial hyperbolic orbit, and where the 
final equality defines x = b/bp. The criterion for obtaining a bound orbit (from the previous section) is equivalent 


3 Note that the sphere of influence, as defined here, corresponds to the location where the incoming trajectory switches from a two-body 
problem with central mass M to a two-body problem with central mass m in the matched conics approximation. The boundary Rgor is 
comparable to, but not equivalent to, the Hill radius Ry = a(m/3M)?/3. 
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to the requirement b < by (x < 1). Since the cross section depends on b?, the distribution of impact parameters will 
be weighted towards larger values. This finding, in turn, implies that typical bound orbits will have final semi-major 
axes comparable to the starting (negative, hyperbolic) semi-major axis of the incoming orbit. For væ = 1 km/s, for 
example, bound orbits will typically have o ~ 1000 au. In order to obtain tighter orbits comparable to the size of the 
solar system (or even the Kuiper belt), we need ay < 100 au, which in turn implies that b bo /10 ~ 10 au. 

If we assume that the impact parameters b are uniformly distributed over an area, with a maximum value bo, then 
the probability distribution for the dimensionless quantity x has the simple form dP = 2ydy. Using equation (26), we 
can determine the probability distribution for the semi-major axes of the bound orbits, i.e., 

dP 2apao 


day (ay + ao)? ` an) 


As written, this distribution is normalized over the interval 0 < ay < Ge, 

Note that the distribution of equation (27) corresponds to the semi-major axes of the bodies when they are captured. 
The orbital elements of the captured objects will continue to evolve (e.g., through continued close encounters with the 
planets), so that quasi-stable orbits will display a different distribution (which should be explored further in future 
work). 


3. NUMERICAL RESULTS 


The cross sections derived in the previous section made use of a number of approximations. In this section we use 
a suite of more than 500 million simulations to numerically compute the capture cross section. 


3.1. Simulation Details 


We sample rocks of mass 107° Mo isotropically on the sphere at a barycentric distance of 10° au. Each rock’s velocity 
unit vector is uniquely defined by its position on the sphere, pointing directly toward the solar system barycenter. We 
then assign each rock an impact parameter at some random angle in its plane tangent to the sphere. We randomly 
sample the impact parameter uniformly given the condition that the maximum pericenter distance qmax < 12 au— 
comfortably above the largest pericenter distance for capture not attributable to chance close encounters with a giant 
planet. Finally we scale the rock’s velocity unit vector by a factor 


2 
ue da, + E (28) 
r 


where u = G} ;m; and i € {Sun, Jupiter, Saturn, Uranus, Neptune}. In Equation 28, vx is the rock’s field (or 
cluster) velocity at infinity, and the second term accounts for the kinetic energy that the rock gains by falling from 
infinity to a barycentric distance r. 

The above procedure gives us a state vector, from which we compute a body’s Keplerian orbital elements. To save 
computation time, we use these elements to propagate each rock along its unperturbed hyperbolic orbit to a barycentric 
distance of 1,000 au. This approximation (that the solar system is a point mass with all of its mass at the barycenter) 
should be accurate to about one part in 10°, since the solar system’s quadrupole term goes like r73. Once we have 
performed the analytic propagation of the rock, we use NASA’s development ephemerides to initialize the solar system 
at a random date in a 200-year range around the arbitrarily chosen Julian Date 2459010.5. This ensures that our 
results are not affected by some exceptional coincidence in the initial phases of the giant planets’ orbits. 

When we have initialized our rock and the solar system, we use Rebound’s IAS15 integrator (Rein & Spiegel 2015) 
to evolve the system numerically. For each simulation, we conserve the system’s total energy to better than one part 
in 10'4—much smaller than the fraction of the system’s energy attributable to the rock. Therefore we are confident 
in the accuracy of our integrations. 

For each integration, there are three possible outcomes: the rock may be captured; undergo a collision with another 
body; or be ejected from the system. If at any point during the simulation the rock’s energy drops below zero, we 
consider it to be captured and end the simulation. If the rock undergoes a collision or if the rock is unbound and 
exiting the solar system with a barycentric distance greater than 40 au, we end the simulation and determine that the 
rock was not captured. We then follow up on our captured objects, integrating for 51% of an orbital period to ensure 
that each object is truly bound (as opposed to having a transient bound osculating semi-major axis due to the phases 
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of the giant planets). If during our followup the object’s apocenter distance exceeds 1 parsec, we consider it to be lost 
to cluster or galactic tides. 

Current models of solar system formation predict that the giant planets formed in a more compact arrangement, and 
then migrated to their current orbits. To account for this we ran a set of simulations with the compact solar system 
model presented in Tsiganis et al. (2005). The cross section we calculate with this model differs from that calculated 
using the solar system at the current epoch by less than 1 percent, so our calculations should be equally applicable to 
the pre-and-post-instability architectures of the solar system. 


3.2. Capture Cross Section 


Since we sampled events uniformly in impact parameter, we can calculate the capture cross section as 
Zoch 
o = 2 biði (29) 


where bmax is the maximum impact parameter sampled, N is the number of events, and 6; is a Kronecker delta that is 
1 if the event resulted in capture, and 0 otherwise. We display our results in Figure 2. As we expect, elt goes like 
v2 in the low-speed limit, and like v% in the high-speed limit. To facilitate the use of the cross section in analytic 
calculations we fit o(v..) with the simple function 


Ou 


Per FI ie 


O(Voo) = 
where u = Uso d'Us and vo is a velocity scale determined by the properties of the planet ejecting the rock. We find the 
data are best fit by parameter values ou = 232, 250 au? and vo = 0.4179 km/s. Keep in mind that these cross sections 
apply for capture into any bound orbit. 

The scale og for the cross section obtained from fitting our numerical results can be compared to the analytic 
estimates of the previous section. If we evaluate equation (19) in the high speed limit, then agreement between the 
analytic and numerical estimates implies that oo = 647a7a7(m/M)?(v,/uz)®, where vy is the orbital speed of Jupiter. 
The expressions are equal if the dimensionless parameter o ~ 0.21. The analytic and numerical results are in agreement 
for all incoming speeds if we identify the scales v, and vo, which is equivalent to considering captures with a maximum 
(post-encounter) semimajor axis dmax * 5090 au. Notice also that vg ~ vz ~ (Gm/a;)\/? (see also Appendix C). 
Similarly, equation (22) agrees with the numerical result in the high speed limit if the dimensionless parameter y ~ 5.8. 

Note that the ensemble of numerical simulations is confined to speeds v3, < 15 km/s. This upper limit is invoked for 
a number of reasons. Due to the steep power-law fall-off of the capture cross section, relatively few capture events take 
place at higher speeds, so additional computation leads to diminishing returns. In addition to the steep dependence 
with va, the numerical data indicate that the power-law begins to break at a comparable speed. Some type of break is 
expected: For væ greater than ~ 10 km/s, capture by close encounters with the Sun becomes ineffective (see Appendix 
A). It is noteworthy that the capture cross section at Us. œ% 13 km/s is comparable to the geometrical area of the Sun 
(~ 7 x 107 au’). For larger encounter speeds, incoming rocky bodies are thus more likely to collide with the Sun 
than be captured into a bound orbit.* 


4. ANALYSIS OF CAPTURED OBJECTS 


In this section we examine the orbital elements of our captured objects to gain insight into the mechanics of the 
capture process. In Figure 3, we show the impact parameter (and pericenter distance) distribution of the unperturbed 
orbits of our captured objects for asymptotic speeds vœ of 1 and 2 km/s. Each histogram displays a clear relative 
peak at the pericenter distances corresponding to the orbit of Jupiter, along with a much smaller peak for the orbit 
of Saturn. Comparison of the two histograms indicates that somewhere between 1 and 2 km/s, the dominant capture 
process switches from that due to the motion of the solar system barycenter to close encounters with a giant planet 
(especially Jupiter). 

In Figure 4, we show the post-capture eccentricity as a function of semi-major axis for a subset of the captured 
objects with vœ of 1 km/s (top panel) and 2 km/s (bottom panel). The figure also includes equi-pericenter curves 


4 For completeness, we note that due to gravitational focusing, the collision cross section with the Sun is larger than the capture cross section 
for speeds væ > 2 — 3 km/s. 
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Figure 2. Capture cross section (in au”) as a function of the asymptotic speed væ (in km/s). The black points represent 
the numerically calculated cross sections, and the corresponding error bars represent one standard deviation. The red curve 
represents the best fit of equation (30) to the data. The blue dashed line shows a power-law of the form o ~ v3, as expected in 
the limit of high velocity. The green dashed line shows a power-law of the form o ~ v3’, as expected in the low velocity limit. 
The shaded regions indicate the parameter space where captures due the Sun (green) and planets (blue) dominate, although 
the boundaries are not sharp. 


corresponding to integer multiples of the spheres of influence of Jupiter and Saturn. The numerical results for captures 
display a relative overdensity of points with pericenter distances at Jupiter and Saturn, indicating that these captures 
are (likely) attributable to close encounters. 
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Figure 3. Histograms of the (unperturbed) impact parameter distribution of captured objects for asymptotic speeds væ of 1 
km/s (top) and 2 km/s (bottom). For convenience, we also indicate the pericenter distance of the unperturbed orbit. Both plots 
show relative peaks at pericenter distances corresponding to the orbits of Jupiter and Saturn. As væ increases, close encounters 
with the giant planets become more important for capture. 
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Figure 4. Post-capture eccentricity versus semi-major axis of captured objects for incoming speeds of 1 km/s (top) and at 
2 km/s (bottom). The orange and red regions correspond to integer multiples of the radius of influence centered at the the 
semimajor axis of the orbits of Jupiter and Saturn, respectively. 


In Figure 5, we show the kernel density representations for the post-capture inclination semi-major axis, and eccen- 
tricity for væ = 0.5, 1, and 2 km/s. While captures become increasingly rare for higher-velocity events, the resulting 
semi-major axes of the captured objects are typically smaller than those for objects captured in low-velocity events. 
This trend is important for assessing object retention, as captured bodies with semi-major axes a = 1000 au are more 
likely to be stripped from the solar system by interacting with passing stars (in the solar birth cluster) or the galactic 
tides (in the field). 
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It is noteworthy that capture events readily produce highly-inclined and even retrograde objects. This finding 
indicates that capture is yet another potential channel for the production of the observed populations of highly- 
inclined and retrograde centaurs, which are currently best explained by the putative Planet Nine (Batygin et al. 2019). 
However, it is important to note that the orbits of the captured objects will evolve over time. As a result, the captured 
objects do not represent a long-term stable population. As the orbits of the captured objects evolve, some will become 
more eccentric until they collide with the Sun; some will undergo scattering events or interactions with the Galactic 
tides and be ejected from the solar system; some will be frozen into the inner Oort cloud by passing stars; and some 
will continue to evolve on stable or quasi-stable orbits (either by becoming caught in resonances with the giant planets, 
or by achieving orbits that otherwise avoid close encounters). 

Because most capture events resulted in high-eccentricity orbits, we have rather low statistics for small values of 
eccentricity. Despite the data limitations, though, it is clear that as vs increases, the low-eccentricity tail of the 
distribution becomes fatter. 
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Figure 5. (Top) Gaussian kernel density estimate of the post-capture inclination distribution (with 7 measured from the ecliptic 
plane) of captured objects at speeds væ = 0.5, 1, and 2 km/s. (Center) Gaussian kernel density estimate of the post-capture 
semi-major axis distribution of captured objects for væ = 0.5, 1, and 2 km/s. (Bottom) Relative fraction (in logarithmic scale) 
of the post-capture eccentricity distribution of captured objects for væ = 0.5, 1, and 2 km/s. Note that these curves represent 
probability distribution functions; there will be fewer captures in the high-speed case than in the low-speed case. 


5. APPLICATIONS 
5.1. Velocity Averaged Cross Sections 


This paper determines the velocity dependent cross section o (Və), which can be fit with a function of the form 


= where us, (31) 


alvo) = Tar ae Se 
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The capture rate for rocky bodies by our solar system is given by 
rT = nrl(ovo) ; (32) 


where np is the number density of rocks that the solar system encounters. The capture rate depends on the velocity- 
averaged cross section, which is given by the integral 


i= f E AE ibs: (33) 


where f(v) is the distribution of encounter velocities of the rocky bodies. 

The distribution of relative speeds f(v) depends on the environment. In the solar birth cluster, f (və) is determined 
by the processes that eject the rocky bodies from their original planetary systems. In general, the clusters are not 
sufficiently long-lived for the rocks to attain a thermal distribution of speeds. Instead, they are expected to retain the 
velocity distribution resulting from the ejection process. If ejection occurs through scattering interactions with giant 
planets, then f(v..) takes the approximate form 


Lët = 


4U /Up 
E /v2) : 
where the velocity scale v2 ~ GM, dos, where ap is the semi-major axis of the planet that scatters the rocks (e.g., see 
Moorhead & Adams 2005 for a derivation). 

Note that the distribution (34) is normalized over the entire interval 0 < væ < co. In practice, the distribution will 
have a maximum value determined by the escape speed from the planets that scatter the rocky bodies. Notice also 
that the full distribution will be a convolution of the distribution of ejection speeds from each planet that scatters 
rocky material. As an approximation, we consider only a single distribution and interpret the velocity scale vp as a 
typical value. As a result, vp is expected to be comparable to the orbit speed of outer planets, i.e., vp ~ 10 km/s. 
Finally, we are assuming that equation (34) corresponds to the distribution of relative speeds between the rocks and 
the solar system (e.g., see the discussion of Binney & Tremaine 2008). 

Putting the above considerations together, we can write the velocity averaged cross section in the form 


e PS Adu v3 
E J e a: (35) 
0 Vv. 


el, UPEP Pe 


(34) 


where u = oe d'Us (as before), we have defined 7 = Ge /Vp, and where the second equality defines the integral function 
I(n). The dimensionless function I(7) can be evaluated to obtain 


3 
oe (1+ TN 37° 44 ) 


Note that I — ~v in the limit 7 — 0, and in practice 7 ~ 1/10. As a result, a good approximation for the capture cross 
section takes the form 


(36) 


a 
(Tv0) & Ou. (37) 

D 

P 
For comparison, we can determine the velocity averaged cross section for the scenario where the cluster rocks are 


virialized and have the same (Maxwellian) velocity distribution as the stars.” In this limit, (ovæ) can be written in 


the form 
Al exp[—v2, Zei 
EE (38) 


where s is the velocity dispersion of the distribution. Note that the value of s for the distribution of relative speeds 
is larger than the value so for the velocity distribution of the stars in the clusters (s = v2sọ). Here we define the 
variable w = Ga dë and the parameter € = s/v,, so that 


lova) = 0022 [wa Pe =o 216), (39) 


5 In general, we expect the stars to reach virial equilibrium much faster than the rocky ejecta. The stars start their cluster trajectories with 
sub-virial speeds, but then fall toward the cluster core where interactions take place, and equilibrium is rapidly realized (in a few Myr, e.g., 
Adams et al. 2006). In contrast, the rocks are ejected with speeds much larger than the virial speed and have little chance for interactions to 
slow them down. Moreover, the stellar virialization starts as soon as stars form, whereas the planet formation and the subsequent ejection 
of rocks occurs many Myr later. 
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where the second equality defines the integral function J(€). The exact form for J(€) can be found. If we define 
u = 1/(2€2) = v2 /(2s?), then 
JE = J(u) = wll - ue" Ei (uw) , (40) 


where E(x) is the exponential integral (Abramowitz & Stegun 1972). In the limit € > 0 (u — oo), the function 
J = 1; in the opposite limit € >> 1 (u > 0), J = u = 1/(2€7). A good working approximation for the cross section of 


equation (39) thus has the form 
4/2 
Vo Vin (41) 


which is exact in the limits and has a relative error less than ~ 20% over the entire range 0 < € < oo. One can use 
equations (39) and (40) if higher accuracy is required. 


5.2. Rock Capture in the Birth Cluster 


Using the results derived above, we can estimate the total mass in rocky bodies that were captured while the Sun 
remained in its birth cluster. The capture rate is given by equation (32), with velocity-averaged cross section specified 
through equations (39) and (40). The rocky bodies will have a distribution of sizes g(R), which is defined here such 
that 


mp = [ g(R)dR and PR= LR g(R)m(R)dR, (42) 


where m(R) is the mass of the rock as a function of its size. With these definitions, the capture rate I can be converted 
into a mass accretion rate given by 
M= PR(OVx~) b (43) 


where pr is the mass density of the cluster in the form of rocks. Given that each planetary system in the cluster is 
expected to eject a few Earth masses of rocky material (e.g., Rice & Laughlin 2019), the density pr is given by 
aMsN. 
pr = —°— =aMon., (44) 
V 

where o is a dimensionless factor of order unity and n, is the number density of stars. For completeness, note that the 
inclusion of icy planetesimals will increase this density estimate. In any case, the total mass in rocky bodies captured 
by the solar system during its cluster phase can be written in the form 


(AM)r = als | | > neat Lem) = AMg (nx) T(oVs0) - (45) 


The final equality defines the mean density of the cluster, where 7 is its effective lifetime. A number of studies 
have found upper bounds on the product (n,)7 by requiring that the solar system is not overly disrupted, including 
considerations of the planetary orbits (Adams & Laughlin 2001; Adams 2010; Li & Adams 2015), the Kuiper Belt 
(Moore et al. 2020), and the orientation of the plane of the cold classicals (Batygin et al. 2020). This work indicates 
that the product is bounded by (n,)7 < 2 x 10* pe? Myr. 

If we take vg = 0.5 km/s, vp = 10 km/s, and ou = 2 x 10° au’, then the velocity averaged capture cross section 
becomes (av..) zs 800 au? km/s. With this cross section, the total mass in captured rocky material from equation 
(45) is about (AM)r ~ LU "Ms, Of course, most of this material will be ejected back into the cluster or the field. 
The retention rate of material in the inner Oort cloud is ~ 1% (Brasser et al. 2006), so we would expect ~ 10-°Mg 
to be captured in the inner Oort cloud. 

Note that if rocks ejected from planetary systems in clusters follow the velocity distribution of equation (34), then 
some fraction of the material will leave the cluster during its first crossing. The high speed tail of the velocity 
distribution will thus be de-populated. In practice, however, most of the capture events arise from the low-speed 
portion of the distribution, so that the correction for the loss of high speed material is modest. 

We can also estimate the mass of rocks captured while the solar system is in the field. In this case, we expect the 
rocky material to encounter the solar system with a velocity distribution comparable to that of the field stars, i.e., 
a Gaussian distribution with s ~ 40 km/s. In this case, the velocity averaged cross section Los) © 0.08 au? km/s. 
If we also assume that each planetary system ejected the same mass in rocks during its formative phases, then the 
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density of rocky material will be proportional to the stellar density (we are thus assuming negligible losses). As a 
result, the product (n.)T ~ 460 pe~? Myr, and the expected mass in captured rocks is about (AM)p © 2 x 107° Mọ. 
Using the approximate retention rate of 1%, we would expect only ~ 2 x 1071! Mg of these rocks to remain in the 
inner Oort cloud. This inventory of captured alien material from the field is exceedingly small, roughly the equivalent 
of one 5 km body. Rock capture during the birth cluster phase is thus expected to produce the dominant contribution 
(by roughly a factor of one million). These latter objects are expected to have radiogenic ages comparable to ordinary 
solar system bodies, but might be identified by different (unusual) chemical composition. 

Note that the values presented here are highly approximate. Not all of the rocks will be captured in the inner Oort 
cloud, so that the retention fraction could be smaller than assumed here (most of the captured interstellar bodies 
initially have Jupiter-crossing orbits, whereas the planetesimals in the Oort cloud could have different origins). In 
any case, most of the captured objects will be ejected, and some will eventually collide with the Sun. Although these 
calculations provide working order-of-magnitude estimates, in forthcoming work we will refine these projections by 
numerically investigating the long-term behavior of the captured bodies from this work. 


6. CONCLUSIONS 


This paper has revisited the problem of capturing interstellar objects on initially hyperbolic trajectories into bound 
states. Using an ensemble of 500 million numerical fly-by simulations, the main result of this study is the determination 
of our solar system’s capture cross section as a function of encounter speed (see Figure 2). The resulting capture cross 
section shows the power-law velocity dependence o ~ v;? in the limit of low speeds and the dependence o ~ v° in 
the limit of high speeds. The capture cross section o(vg,) over the entire range of asymptotic speeds can be fit with 
the function given in equation (31). 

This paper also presents an analytic treatment of the capture problem using the approximation of matched conics and 
the (inverse) gravitational slingshot effect (Section 2). These arguments show that capture by both close encounters 
with the Sun (Section 2.2), and by close encounters with a giant planet (Section 2.3), have the same nearly velocity 
dependence as that seen in the numerical simulations (namely o ~ v}? at low speeds and o ~ v8 at high speeds). 

The capture events can be classified as either close encounters with the Sun or close encounters with giant planets. 
At low speeds, encounters with the Sun dominate the capture cross section. At higher speeds, close encounters with 
Jupiter dominate. Close encounters with the other giant planets contribute to the cross section, but do not dominate 
the dynamics. More specifically, for the particular case of veo = 1 km/s, capture events due to close encounters with 
Jupiter are ~ 100 times more likely than captures due to Saturn. The frequency of close encounters with Uranus and 
Neptune are smaller (than for Saturn) by an additional factor of ~ 100. 

With the capture cross section as a function of velocity specified, the effective mean cross section (va)/(v) can 
be determined for any distribution of encounter speeds. For the case of a Maxwellian distribution and a power-law 
distribution motivated by rock ejection, the mean cross section can be evaluated analytically (see Section 5.1). 

Finally, as an application of the capture cross section, we estimate the total mass (AM), in the Oort cloud that 
originates from other planetary solar systems (Section 5.2). The mass accreted while the Sun lived within its birth 
cluster is of order (AM)p ~ 107° Mẹ, about a million time larger than the mass subsequently accreted from the field. 

Although the capture cross section for the solar system is now well-characterized, many avenues for future research 
remain. The simulations of this paper consider the capture of interstellar objects and the resulting cross sections 
include all capture events, independent of their residence time in the solar system as bound objects. Future work 
should determine how long captured bodies can remain bound to the Sun, since many such objects are expected to 
be ejected from the system or to collide with other solar system members. The residence time (ejection time) should 
thus be determined for each type of orbit displayed by the captured objects. With these results in place, one can make 
a refined estimate of the current population of alien objects in the solar system, along with their expected orbital 
properties. 
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APPENDIX 


A. UPPER BOUND ON INCOMING SPEED FOR CAPTURE 


In this Appendix, we find upper limits on the asymptotic speed, i.e., the largest speed for which a rock can be 
captured by the solar system. We consider both close encounters with the Sun and with one of the giant planets. 
We start with the capture criterion of equation (13), which we reproduce here, 


v2, < 4Uvcos@ — 4U? < An, (Al) 


where the second inequality follows because cos @ < 1 and 4U? > 0. The largest speed U possible for the Sun, relative 
to the center of mass, is the orbit speed due to the giant planets, where 


1/2 1/2 
yam G(M +m) ~ D GM . (A2) 
M ap M \ ap 


where the second expression is consistent with the ordering approximation and where we have included the reflex speed 
due to only one planet. The maximum possible speed is given by the sum of the contributions of all of the planets. 
We can account for this complication by writing the limit in the form 


Zen: COMM? 
GE E l KS 
Thus far, our upper bound has the form 
1/2 1/2 1/2 
2 mJ GM — MJ GM 2 2GM 
UT Bur (=Z ) v= Bur (Z Uso + Ro ; (A4) 


where we have taken the minimum distance (and maximum speed) of the orbit to be given by the radius of the Sun. 
If we drop the first term in square brackets, the expression simplifies to the form 


GM\!P 2G6GM\ 1? Gm 
SE CR A5 
Us, < ver Ro VIRB (8km/s) (A5) 


To derive a rigorous version of the upper bound, we define two velocity components 


_ omy (GM\'? _ (2GM\*? 
UE Sr E and w= Ro e (A6) 
so that the limit takes the form 
S vi 1/2 ue 
Us, < uw +o Eog? (A7) 


Since u is of order the orbital speed of the Sun (specifically, a small fraction of 1 km/s) and w is the escape speed of 
the Sun (~ 620 km/s), we find u < w. In this limit, the full expression of equation (A7) reduces to the form v2, < uw, 
which corresponds to the approximation of equation (A5). 

We can also consider the case where the incoming rock enters into the sphere of influence of a planet and loses energy 
through an inverse gravitational assist from the encounter. In order for the rock to be captured by the solar system, 
the final speed must be sufficiently small, i.e., 

Ue < ou (A8) 
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where the orbital speed vor» of the planet is a measure of the depth of the gravitational potential well at the location 
of the planet (the location of the encounter). We also assume that the incoming rock obeys conservation of energy so 
et 1 1 GM 1 
2 2 a} 2 2 2 
ie = E = 3” „3Y T Vorb => UW = Dës + 2U5 +5 - (A9) 


Using this expression in the result for the post-encounter speed from before, we obtain 


v% + AUT < 4U cos 6[v2, + 2v2]. (A10) 
If we take U = vor» cos ¢, this expression can be rewritten in the form 
vs, < 8u24,02, cos” 6(2 cos” 0 — 1) + 164, cos? 6(2 cos” 0 — cos? dl. (A11) 


We can thus obtain an upper limit on the asymptotic speed for an object to be captured by taking cos 0 = 1 = cos @, 
i.e., 
vt, Baatz, + 16v4,,, (A12) 


which leads to the bound 
Væ < Worp(1 + V2)? (A13) 


The fastest orbit is that of Jupiter, where Vorb ~ 13 km/s, so we have the limit va < 40 km/s. 

We thus find that the maximum possible speed (for capture) is larger for the channel involving close approaches 
to planets (in particular Jupiter) than for close approaches to the Sun. On the other hand, planet encounters are 
expected to occur much less frequently. Note that these limits do not account for effects such as radiation pressure 
or atmospheric drag. It is therefore possible in principle to capture rocks with larger væ, but such events would be 
exceptionally rare. 


B. ROCK CAPTURE BY CIRCUMSTELLAR DISKS 


Rock capture could also take place due to gas drag if interstellar bodies enter the solar nebula while it retains its 
gaseous component. This Appendix explores the efficacy of this process. 
If a rock has speed vo when it enters the gaseous region of the disk, it will have final speed ur given by 


vf = vo exp|—pAl/m] , (B14) 


where p is the gas density, A is the rock area, m is the rock mass, and £ is the length traveled. If a rock comes from 
infinity with speed væ, then it will have a larger velocity when it hits the disk, so that vo is given by 
2GM 
2 EE 


vg =v (B15) 


2 r 


In order for the rock to lose enough energy to enter into a bound orbit, the final speed must be less than the limit 


2GM 2pAl 2GM 
v= D if Jew | p < . (B16) 
r m r 
The capture criterion thus becomes 
2pAl Gr 
ECH 1 Sea P B17 
SS EE | ý Kg Sg 
If the rock passes through the disk vertically, then of = X, where the surface density for the solar nebula can be 
written in the form 
lau\?? 
U(r) = än (=) (B18) 
r 


where ©, œ 3000 g/cm?. If the rock passes through the nebula at an arbitrary angle 0, then pf = X/cos@. We can 
also write A = mR? and m = pr(47/3)R? where R is the size of the rock. If we work in the low speed limit with væ 
= 1 km/s, and write R in units km and r in units of au, then the criterion for capture reduces to the simple form 


> Bräi. (B19) 
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For the fiducial case where cos 0 = 1, the solar nebula can thus capture incoming rocks with radii R <1 km at radial 
location r = 1 au. Somewhat larger rocks can be captured for typical inclined trajectories (e.g., the limit becomes 
R<2 km for cos = 1/2). Smaller rocks are readily captured, and rocks as large as dwarf planets are highly unlikely 
to be captured (see also Brasser et al. 2007). For incoming trajectories that are confined to the plane of the disk, the 
above treatment must be modified to include the disk structure (i.e., the above treatment does not apply in the limit 
vos — 0). 


C. DIMENSIONAL ANALYSIS 


The capture problem has more than one dimensionless field, so that the cross section of interest cannot be directly 
determined from dimensional analysis. Nonetheless, such an analysis is useful for understanding the result. 

Consider the solar system to consist of only the Sun-Jupiter binary. The physical variables required to characterize 
the systems are thus the masses (M, m) and the orbital radius of Jupiter az. Note that the eccentricity of the Jovian 
orbit is too small to matter, and the radii of the bodies also do not contribute. The other variables in the problem are 
the gravitational constant G and the speed of the incoming rock o, Here we want to determine the cross section as 
a function of these variables. We thus want to determine the impact parameter b of the incoming orbits that allow for 
capture. 

The first estimate for the cross section is thus 

oy = 707. (C20) 


The remaining (relevant) variables can be combined to construct two dimensionless fields, i.e., 


A = = and A2 = an (C21) 


2 2 i 
Vso AI Vso aI 


Alternatively, one could define the second field to be the mass ratio m/M, and the quantity Ag = Ay(m/M). At low 
speeds, the capture cross section has the form 
02 = mas Ay ; (C22) 


whereas at high speeds the cross section becomes 
o3 = taf A A2. (C23) 


We can combine all of the above results to write the capture cross section in terms of the dimensionless fields of the 


problem: 
2 MA3 


04 ~ BEIT Aa)? | (C24) 
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